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Previously published CTEQ6 parton distributions adopt the conventional zero-mass parton scheme; 
these sets are most appropriate for use with massless hard-scattering matrix elements commonly found in 
most physics applications. For precision observables which are sensitive to charm and bottom quark mass 
effects, we provide in this paper an additional CTEQ6HQ parton distribution set determined in the more 
general variable flavor number scheme which incorporates heavy flavor mass effects. The results are obtained 
by combining these parton distributions with consistently matched DIS structure functions computed in the 
same scheme. We describe the analysis procedure, examine the predominant features of the new distributions, 
and compare with previous distributions. 
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I. INTRODUCTION 

Parton distributions provide the essential link between the theoretically calculated partonic cross- 
sections, and the experimentally measured physical cross-sections involving hadrons and mesons. This link 
is crucial if we are to make incisive tests of the standard model, and search for subtle deviations which might 
signal new physics. Since perturbative calculations are renormalization scheme dependent, "'^ it is important to 
use properly matched hard-scattering cross-sections and parton distribution functions (PDFs) in evaluating 
factorized cross-sections for physical applications. This issue is particularly relevant for applications involving 
heavy quarks, since the heavy quark introduces a new mass scale which leads to complications of the PQCD 
formalism. 

For most physical applications, it is convenient (and a good approximation) to use hard-scattering 
cross-sections calculated in the zero parton-mass limit. In fact, next-to- leading-order hard-scattering cross- 
sections in the non-zero quark-mass case have only been calculated for some very basic processes, such as 
deep inelastic scattering (DIS); they are not yet available for most interesting "new physics" applications. 
Thus, the most useful parton distributions for general applications are those determined in a global QCD 
analysis using the zero quark-mass approximation. This was the choice made for the latest series of CTEQ6 
parton distributions: CTEQ6M, CTEQ6D, and CTEQ6L. [3| For definiteness, this scheme will be referred 
to as the zero-mass variable-flavor- number scheme (ZM-VFNS). 

This paper extends Ref. by performing a similarglobal analysis using the generalized (non-zero 
quark-mass) MS perturbative QCD framework of Refs. 0, IE ^ IE Si which we label the general-mass 
variable-flavor- number scheme (GM-VFNS). When matched to the corresponding hard-scattering cross- 
sections calculated in the same scheme, the combination should provide a more accurate description of 
the precision DIS structure function data, as well as other processes which are sensitive to charm and bot- 
tom mass effects. The main result we present here is a new set of parton distributions evaluated in this 
GM-VFNS which we identify as CTEQ6HQ.2 

Perturbative QCD with non-zero quark-masses is well-established in the literature. 13. Isll^ Eolllll 
ll^ Its implementation is, however, necessarily more complicated than the corresponding zero-mass case. In 
addition to the familiar renormalization-scheme dependence, more ambiguities arise from the implementation 
of the charm and bottom thresholds (in relation to their masses), and the detailed way the 3-flavor, 4-flavor, 
and 5-flavor renormalization schemes are matched at the transition scales. Several schemes proposed in 
the recent literature are formally equivalent to each other at high energies; but, they can differ in the 
threshold reg ion, as well as in the intermediate energy range, where much of the experimental data lie. 
Qi IE llll- I12I I13I Il4j | Some of these schemes are more natural and numerically robust, e.g. insensitive to 
variations of the remaining parameters (such as scale), than others. We follow the procedure proposed in 
Refs. [313 which has been shown to be particularly stable as it naturally takes into account the dominant 
(logarithmic and non-logarithmic) mass-dependent kinematic effects with a physically motivated choice of 
the scaling variable. 

Section 2 summarizes the method of analysis. Because this work represents an extension of the zero 
quark-mass CTEQ6 global QCD analysis, we shall focus mostly on the new elements related to non-zero 
quark mass; details that are in common with the massless case can be found in Ref. 0. Section 3 describes 
the global analysis of data. Section 4 compares the new parton distributions with data. Section 5 makes 
comparisons between the CTEQ5M, CTEQ6M and CTEQ6IIQ PDFs. At asymptotically high energies, 
mass effects become negligible; therefore, all generalized parton model prescriptions should approach the 



^ We shall use the term renormalization scheme in the general sense to include both the usual renormalization (ultra-violet 

subtraction) scheme, and the factorization (infrared and coUinear subtraction) scheme. 
2 A fit named CTEQ6F3 employing the three fixed-flavor-number scheme (3-FFNS) for DIS (cf., Refs. |3|3|) which is suitable 

for special applications in that scheme, will be presented separately. 
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same zero quark-mass limit. We show that the differences between CTEQ6M and CTEQ6HQ PDFs indeed 
vanish with increasing Q, as expected. This section also contains a discussion of applications and issues 
related to uncertainties of physics predictions resulting from these parton distributions. 



II. THE GENERALIZED MS SCHEME WITH NON-ZERO MASS HEAVY FLAVOR PARTONS 

The main feature of this analysis, compared to the standard zero-mass parton CTEQ6 PDF sets, is the 
adoption of the generalized perturbative QCD framework incorporating non-zero mass effects for the charm 
and bottom quarks.'^ Collins has shown that, to all orders of the perturbation expansion, factorization of 
the DIS structure functions holds for massive partons to the same degree of rigor as in the familiar zero- 
mass cas e, la The relatively simple physical ideas underlying the full factorization proof are described in 
Refs. [il Il4i lisf . To implement this formalism, which we shall refer to as the generalized MS scheme (or 
GM-VFNS, as mentioned earher), we follow the exphcit procedures outlined in the Appendix of Ref. 0], and 
supplement this with the threshold prescription of Refs. 0, . The latter reference also contains a review 
of alternate approaches. 

The NLO DIS structure functions are given by the generic formula (suppressing the structure function 
label 1,2,3), 

n^-W^E/f /.^^""(f.f.^.M/.)) (1) 

where "a" is the initial state parton label, fa is the parton distribution function, lu"" is the hard-scattering 
cross-section calculated in PQCD to order al, A is the QCD-lambda parameter, and niH (generically) 
represents heavy quark masses, if present. The prescription-dependence allowed by the PQCD formalism 
is associated with possible implementations of the first term on the right-hand side of Eq. ^ within the 
accuracy specified by the remainder term (which will be dropped from now on)."* In the sum over partons in 
Eq. ^ contributions due to the light quark (u, d, s) are standard; contributions arising from massive charm 
and bottom quarks, along with that of the gluon, are treated differently, as explained below. 



A. Leading contributions and the ACOT(x) implementation 

To be specific, we will consider DIS neutral-current charm production process. The leading-order 
process is a virtual-photon scattering off a charm parton: 7*c c. The NLO QCD corrections for this 
process are the boson-gluon fusion process (7*5 — > cc), and the gluon radiation process (7*c — > eg) as well 
as the charm- initated one- loop virtual process (7*0 — > c). 

First, let us examine the kinematic region near the charm production threshold: W = ^Q'^iXjx — 1) ~ 
2 nic- To make the notation precise, we will use the renormalization scale /i to generically label our charac- 
teristic energy scale; separately, we use [im to label the matching 'point, and /ij to label the transition point 
(cf. ^ and below). To implement the generalized MS scheme, the first step is to consider a 3-fiavor scheme 
(appropriate for energy scales fj,^mc), a 4-fiavor scheme (for scales fi^nic), and choose a matching point 
firn (on the order of rric) where the two schemes are matched (i.e., where the discontinuities of as{fi) and 
fa{x,pL) are calculated). We choose the conventional value for the matching point, — ^c, which yields 
the simplest matching conditions. 0] 



^ The top quark can, for all practical purposes, be treated as a heavy particle, not a parton. 

* Note, according to the factorization proof of Ref. the remainder term contains corrections of order A? /Q"^ and A? /m?^ 
but no corrections of order m?^/Q^ . 
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Having precisely defined the 3-fiavor and 4-flavor schemes, which co-exist in the region of the charm 
threshold, we can still choose a transition point, /ij, where one makes the transition from the 3-flavor to the 
4-flavor scheme in the calculation of physical quantities.^ The choice of the transition point is arbitrary — its 
choice is part of the definition of the composite renormalization scheme. While there is much flexibility in 
the choice of fit, clearly it must lie within the overlapping region of applicability of the 3-flavor and 4-flavor 
schemes — therefore, we typically choose fit ^ rric. Although a case can be made, in principle, for choosing 
fit far above TOc (since the 4-flavor scheme is not a natural scheme just above rric), or for considering an x- 
dependent transition scale, fitix, mc) (because the threshold W > 2 rric depends on x as well as Q), we make 
the plain choice fit = to simplify the calculation. This choice is reasonable only within a prescription 
for handling x-dependent threshold effects of Eq. ^ which naturally suppresses the 4-flavor contribution to 
the physical structure function when fi rUc and/or W ^2 mc- The specific renormalization prescription of 
Ref. [iJI has precisely this feature, as we shall outline below. 

Detailed descriptions of the various NLO contributions to the the generalized MS scheme were dis- 
cussed in Ref. 0. While we refrain from reproducing them here, we do need to specify the precise way the 
threshold effects are implemented when heavy quarks are involved. Consider first the numerically significant 
NLO boson-gluon fusion contribution (7*5 cc): 

F^'\l*9^cc:) = / ^<7(-,A^)^J-,^,^) (2) 

J z \z fl fi J 

^ «.(.)[/;^»(..).;(f,^)-..-(£;)/'^.(..)p,^=(9 

where the first term on the right-hand side corresponds to the unsubtracted boson-gluon-fusion diagram 
contribution (with the as factor explicitly taken out),^ and the second term represents the subtraction term 
which renders LUg infra-red safe. The variable x is dictated by the kinematics of the boson-gluon-fusion 
partonic subprocess (7*3 cc) to be a generalized scaling (or "rescaling") variable 

X = .(l + ^) (3) 

where x is the conventional Bjorken variable. In contrast, the other scaling variable, C, in Eq. |21can, in 
principle, be chosen arbitrarily provided the following two constraints arc imposed. 

1. C — > a; in the high energy limit. 

2. The same C variable is used here as in the order ("leading-order") process 7*c ^ c. 

In the high energy limit (m^/Q^ — > 0), the first constraint ensures that both —> x and x ~^ ^ such that 
the complete boson-gluon fusion contribution is infra-red safe. The leading-order process 7*c — > c is: 

F(o)(rc-.c)./^c(z,,)^og^:i^) ^(l-^)=c(C,M) . (4) 

For the second constraint, if the same C variable is used for both the boson-gluon fusion process of Eq. Q 
(7*5 ~^ cc) and the leading-order (7*0 c) contribution of Eq. then the gluon subtraction term 



^ The distinction between the matching and the transition points is not addressed in most papers. This distinction is made in 
Ref. it is discussed at length in the Appendix of Ref. (jj. As mentioned above, arguably, there is a good case for choosing 
the transition point fit to be higher than the simple matching point 

® We use the notation that, for each flavor, uja denotes the parton-level cross-section before infrared and collinear subtraction, 
while uta denotes the corresponding infra-red safe cross-section after subtraction. 
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(last term in Eq. [2J) cancels the leading-order contribution in the limit fi rric] this ensures that the 4- 
flavor formula will reduce to the 3-flavor formula, and the dominant contribution in this region will be the 
unsubtracted gluon fusion term, as required by the correct physics to this order of PQCD. 

A natural choice of C, that ensures the best behavior at both the high and the low energy limits is 
proposed by Ref. Q: 

C = x = x{l+Aml/Q^) . (5) 

This choice has the property that as the heavy quark production threshold is approached from above iW — > 
2mc), then C ^ 1 such that the integration measure in Eq.|21 vanishes.^ Therefore, with the choice C = Xi the 
subtraction term in Eq.Eland the LO quark term, Eq.^ individually vanish at threshold by kinematics. In 
addition, the cancellation between these two terms due to the dynamics of PQCD, as explained in Ref. 
still operates. 

This results in unequaled stability of the predictions in the threshold region, compared to previous 
implementations of the GM- VNFS in the threshold region 0, fl^ . These usually give rise to rather 

prescription-sensitive behavior just above threshold because the perturbation-expansion-inspired prescrip- 
tions generally overlook heavy quark production kinematic requirements, while the physically important 
gluon fusion term (in the threshold region) is strongly kinematically-suppressed. The robust implementation 
we adopt will be referred to as the "ACOT(x)" scheme, following Refs. [l3.ll5|.^ 



B. Other contributions to the full NLO calculation 

For the complete order as calculation, we must also include the real and virtual NLO quark-initiated 
contribution to Eq. ^ 

fW(7*c-.c + X) ^ /^c(z,M)^.^f-,^,^) . (6) 
J z \z ^ J 

The hard-scattering cross-section, toKx, has been computed for general masses 01 the ACOT 

scheme For scales around the canonical DIS choice of /i ~ Q, this term is numerically insignificant com- 
pared to the terms considered above. We include it in our calculation to ensure (perturbative) consistency. 

It is also worth mentioning that the calculation of this term, which is computation-intensive in spite 
of the smallness of its numerical contribution, can be considerably streamlined by using a simplifying feature 
of PQCD in the generalized MS scheme: one can set rric — in uil and in Eq. |5J and Eq. Q without 
sacrificing accuracy. This point is discussed in detail in Ref. {18i] where the numerical insensitivity of the 
calculation to the choice of rric is explicitly demonstrated.^ 



C. Remaining scheme and scale choices 

The only remaining arbitrariness of our prescription is now associated with the familiar choice of the 
factorization (and renormalization) scale fj,. The scale dependence of perturbative calculations in general 
arise from imperfect matching between the terms in the truncated perturbation series. Because of the 
near-perfect matching in the ACOT(x) prescription discussed above (in contrast to alternate choices), the 



This is more easily seen from the equivalent formula x = 1 — {I — x){l — Am'^/W^). 
* So named since it is a variant of the original ACOT approach Q, replacing the naive x with the kinematically natural 

rescaling variable x of Eq. |3 

® The simplified scheme [T8i| can also be recovered by investigation of the resummed perturbation series, [ifij 
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dependence of the physical structure functions on the choice of /i is completely negligible (cf., Ref. 14|) 
throughout the energy range of our global analysis. For definiteness, we use /i^ = + m^. In summary, 
the quark initiated terms in the simplified ACOT(x) scheme adopted in this work are: 

1^ d.z -0,1 f X Q \ 

z \z fl fl J 

with representing the standard coefhcient functions for massless quarks in LO and NLO, defined in the 
conventional MS scheme. 

The same treatment is applied to the 6-quark threshold region. Here, the 4-flavor scheme is matched 
onto the 5-flavor scheme for scales larger than rrib, and the production of open bottom particles in the final 
state is added to the inclusive DIS structure functions. In principle, a similar treatment needs to be applied 
for other processes included in the global analysis (currently, Drell-Yan and jet production) using appropriate 
hard-scattering cross-sections calculated in the massive scheme. In practice, this is not done for two reasons: 
(i) NLO calculations of the hard-scattering cross-section involving massive quarks do not yet exist for these 
processes; and (ii) unlike the case of DIS, the experimental errors in these processes are relatively large 
compared to the anticipated differences between the massless and massive scheme calculations. 



- / -ciz,f,)C^ -,- 



III. GLOBAL FITTING AND THE CTEQ6HQ PARTON DISTRIBUTIONS 

The new global fitting performed in the generalized MS formalism follows the same procedure as 
that of the earlier CTEQ6 analysis using the zero-mass parton formalism. The data sets used before are 
supplemented by: (i) the HI F2 ^ set 0| (which was inadvertently left out in Ref. I'D; and (ii) the HI |2o| 
and ZEUS data sets for the structure function with tagged charm particles in the final state. The 
additional HI F2 ^ data set does not have much influence on the new analysis, since one already obtains 
an excellent fit to these data just by comparing them with predictions of CTEQ6M. The F2 data sets are 
quite relevant for this analysis since is sensitive to the charm and gluon distributions, which are tightly 
coupled in the generalized MS formalism.^" 

The parametrization of the non-perturbative parton distribution functions at Qo — "^c = 1-3 GeV 
is the same as in Ref. 0. For this study, we assume as usual that charm partons are entirely "radiatively 
generated" (i.e., through QCD evolution) from the starting scale Qo onward. This assumption is somewhat 
arbitrary, and it is obviously dependent on the choice of Qo- The possibility for having a small component 
of non-perturbative charm 22] at low Q, and its physical consequence will be examined separately. 

As in the previous CTEQ6 analysis, correlated experimental systematic errors are fully incorporated 
whenever available. The best fit obtained with these inputs shall be called CTEQ6HQ, or C6HQ for short. A 
broad measure of the quality of this fit is provided by the overall of 2033 for a total number of 1950 data 
points (xVDOF = 1.04). This is to be compared to a of 1946 for 1811 points (xV^OF = 1.07) in the 
case of CTEQ6M (abbreviated to C6M in the following) |]| . To gain a better feel of how these fits compare, 
we show in Tabled (the 3rd and 4th columns) the overall x^ values and, in parentheses, x^ per data point, 
as well as for the applicable individual data sets. The total number of data points in this head-to-head 
comparison (not including the charm production data points) is 1925. The new C6HQ fit reduces the overall 
X^ by 29 out of ~2000 as compared to the C6M fit. The improvement of this generalized MS result over 
the zero-mass MS result is encouraging, since the generalized MS formalism represents a more accurate 
formulation of PQCD. However, a difference of x^ of 29 is within the current estimated range of uncertainty 



These data sets were not used in the CTEQ6M analysis because F2 is not well-defined theoretically in the the zero-mass 
parton formalism. 
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TABLE L Comparison of the values of the general-mass CTEQ6HQ fit (3rd column) with the zero-mass CTEQ6M 
fit (4th column) . Also included are comparisons with two "mis-matched" cases when GM and ZM parton distribution 
functions are convoluted with the other (wrong) hard-cross-section (5th and 6th columns). The first number of each 
entry is the value, the number in parenthesis is per number of points. Correlated systematic errors, if available, 
are included. 

of PDF analysis. Q, Esl 01 Therefore, the significance of this difference is arguable. We also note that the 
improvement in is spread over most of the data sets: there is no smoking gun for the overall difference. 

The last two columns of Table compare the above results with two possible uses of the PDFs that 
represent a mis-use of PQCD in principle, but occur frequently in the literature in practice, perhaps out of 
necessity. These involve using PDFs obtained in the general-mass scheme convoluted with hard-scattering 
cross-sections (Wilson coefficients) defined in the zero-mass scheme, and vice versa. For the same data 
sets, these mis- matched schemes result in a difference of 420'--^490 in the overall . These are quite large 
differences relative to the tolerances discussed in Refs. 0,E^, and results in clear discrepancies with some 
of the precision DIS data sets, as can be seen in Tabled The lesson is clear: for quantitative applications, 
it is imperative to maintain consistency between the PDFs and the hard- scattering cross-sections. 

IV. COMPARISON WITH DATA 

In Ref . Q , we presented an extensive comparison between the CTEQ6 results and the data sets used in 
the global analysis, including new ways to explicitly account for the correlated experimental systematic errors. 
Since the new C6HQ fit is generally similar to the C6M fit, we shall not duplicate the same comparisons 
for those quantities where the differences are minimal. In subsection IIV Al we compare with the DIS charm 
production data, which were not used in the previous CTEQ6 analysis. In subsection IIV Bl we discuss some 
implications of the neutrino data. 



For ins tan ce, the MRST distributions l25l are obtained in the General-Mass formalism (using the Robert-Thorne implemen- 
tation they are often used in applications convolving with readily available Zero-Mass hard-scattering cross-sections. 
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A. DIS Charm Production Data from HERA 



Because the HI and ZEUS F2 structure function are sensitive to the charm and gluon distributions, 
these win play an important role in extracting the C6HQ PDFs. In Fig. ^ we display comparisons of 
the charm production data from HI and ZEUS [2^ with the C6HQ theory value. For simplicity, we 
combine all the data points from each experiment in a single plot, and we scale the plot by the theoretical 
calculation. On the x-axis, the data points are ordered using x and Q as the primary and secondary sorting 
variables, respectively. The error bars represent the statistical and uncorrclated systematic errors added in 
quadrature. We see that the fits are good with per data point of 0.881 (7/8) for HI, and 0.98 (18/18) for 
ZEUS, respectively. 





FIG. 1: Comparisons of the charm production structure function, F2, with the data from HI (Fig. ^a) (20| and 
ZEUS (Fig. 0b). 26] The j/-axis is (Experiment-Theory) /Theory; along the x-axis, the data points are ordered using 
X and Q as the primary and secondary sorting variables respectively. The error bars represent the statistical and 
uncorrelated systematic errors added in quadrature. 



B. Neutrino Data in Global Parton Analyses 

The NuTeV measurement of the weak mixing angle sin^ (Ref. has recently focused considerable 
attention on the neutrino induced DIS process. 28] In this measurement, the weak mixing angle sin^ 
is extracted from the ratio of neutral current (NC) to charged current (CC) ly and D cross-sections. Even 
before this sin^ 9^ measurement, there have been other long-standing unresolved issues in comparing CC 
structure functions i^2^'^ (measured in 1/ and D scattering) with NC structure functions F^'^ (measured in 
and /i^ scattering) at modestly low-x 10^^). One manifestation of this is that the recently measured 
structure function AxF^^'^^ is not compatible with the QCD predictions. j^^liOl As this structure function 
is particularly sensitive to the heavy quark components, AxF^^'^^ ~ s{x) — c(a;), it is certainly important to 
properly treat the heavy quark mass. While the deviations of AxF!^^'^'^ are at a non-dramatic ^ la level, 
the pattern has been systematic and persistent. 

Therefore, we re-examine the influence of the neutrino structure functions in the current global analysis 
keeping in mind the following issues. First, the high statistics neutrino structure function data of CCFR used 
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here are from the newer "physics model independent" data analysis .31], rather than the earlier ones which 
contained model-dependent corrections. 32] Secondly, meaningful (i.e., quantitative) comparison of NC 
and CC structure functions (which represent different hard-scattering processes) can be made only within a 
common underlying theoretical framework with an accuracy comparable to the experimental uncertainties.^^ 
The proper treatment of charm mass effects, such as those described in the current study, are an important 
part of that theoretical formalism, cf. [s^ . 

First, we ask whether the neutrino structure function, ¥2^'^^ , obtained by the newer model- 
independent analysis is consistent with the muon- induced NC fixed-target structure functions, F2 ^ , and 
the HERA collider structure functions, F| ^, in a consistent NLO analysis. ^"^ Fig. [21shows the comparison 
between the data with the CTEQ6HQ fit (using rric = 1.3 GeV). We see that, even if this data set is 
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FIG. 2: CCFR structure function (from the physics model independent analysis, Ref. j31|') compared to 

CTEQ6HQ fit: (a) low x bins (left panel) show a systematic disagreement; (b) medium to high x bins (right panel) 
show good agreement. 



included in the global fit (cf.. Table P), the measured values in the low- a: bins (left figure) consistently lie 
above the theory curves. The systematic deviation of data points from the global fit is reflected in the 
relatively large x^'s associated with this data set, as seen in Tabled The fit is dominated by the much more 
extensive NC data (from the BCDMS, NMC, HI, and ZEUS experiments), hence one can regard the theory 
curves as representing the NC data properly "corrected" for NLO QCD effects. The pattern of deviation 
of the neutrino data from the average NC predictions at low- a: is similar to earlier comparisons. In other 



For example, direct comparisons of the ratio F^/Fj* to theory (suggested by the "5/18-rule") are meaningful only at the LO 
parto n m odel level, which is not appropriate for the precision data from current experiments. 

Ref. I31I suggests these data are consistent by examining the ratio of the structure functions. However, as explained in 
the text, it is not appropriate to compare the ratio of two different structure functions with theory in QCD beyond LO; a 
complete analysis involving a NLO fit is required. 



11 



words, the discrepancy between CCFR and the other NC measurements at small x persists in the 
context of contemporary global QCD analyses, even when charm mass effects are properly treated as in the 
present analysis. Other possible theoretical sources for this difference (in addition to the NLO contributions 
included in this comparison) have been discussed in Refs. 28, 29, 34], and in the literature quoted therein. 
On the experimental side, preliminary NuTeV measurements jSSj seem to yield improved agreements with 
the NLO QCD predictions at low- a;. It remains to be seen whether the final NuTeV results will resolve this 
problem. 

In our analysis, the strangeness distribution is constrained to be proportional to [u{x) and its 

normalization at the scale Qo is constrained to be 

^ ^ ^ dxx[s{x) +s{x)] ^ 1 
J dxx[u{x) + d{xy\ 2 

as in most current global analysis. In the future, a NLO analysis of the CO charm production data [vs — > 
Cji jjL^^'^X) from CCFR and NuTeV (36| has the potential to determine {s{x),s{xy\ with much more 
precision. Then, a re-assessment of the apparent discrepancy between the measured /S.xF^^^^ and the low- 
lying QCD predictions [2^ '^iH be warranted. An accurate determination of the strangeness of the proton, 
as well as the strangeness asymmetry \s(x) — sjx)] will have important implications for these measurements, 
including the NuTeV anomaly, cf., Refs. |27ll28|. 

Clearly, the interpretation of the neutrino DIS data requires a detailed NLO analysis. Global QCD 
analyses of parton distributions, such as the one performed in this paper, are essential to unraveling the un- 
derlying physics. Progress on both the experimental and the theoretical fronts are needed; this advancement 
is ongoing. 



V. COMPARISON WITH RELATED PDFS: CTEQ6HQ, CTEQ6M, CTEQ5HQ 

The CTEQ6HQ and CTEQ6M fits provide comparable descriptions of the global QCD data in two 
different schemes. Some of the differences in the PDFs arise purely from the choice of scheme. We are partic- 
ularly interested in the differences for the charm distribution, and the closely correlated gluon distribution, 
due to the improved treatment of heavy quark effects in the generalized MS scheme. It is also interesting 
to compare the differences between the earlier (previous generation) CTEQSHQ (C5HQ) distributions with 
the new CTEQ6IIQ distributions; differences between these PDFs are attributable both to new data, and 
to minor differences in the way the theoretical inputs are implemented. We combine these comparisons in 
the figures presented below. 



A. Light partons at the input scale 

Figs. |3i,b show the u-quark and c?-quark distributions from C5HQ, C6M and C6IIQ at the initial 
Q-scale of 1.3 GeV (the charm mass). In order to exhibit the behavior of the PDFs at both large and small x 
values, the x-axis is scaled linearly in -/ir. To give some additional physical insight to these plots, the y-axis 
is chosen to be x^'^/(x, Q), so that the area under each curve represents the momentum fraction carried 
by that distribution. This chosen scaling then focuses on the valence peak, and suppresses the small- a; sea 
excitations. We see that the valence u-quark and d-quark distributions are well constrained by the precision 
DIS data; the scheme dependence is not pronounced. 

Figs. 0^,b show the same PDFs, except the axes are scaled differently in order to highlight the sea 
quark distributions. The differences between the CTEQ5 and CTEQ6 generations of PDFs are a result of 
the improved low-x data. When we compare the HQ and MS type PDFs of the CTEQ6 generation, we 
observe the general pattern that the HQ distributions overshoot the MS ones at low x (~ 10~^), and that 
the distributions meet again at ultra-small x 10~^). We can qualitatively understand the larger sea quark 
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FIG. 3: Comparison of CTEQ5HQ, CTEQ6M, and CTEQ6HQ parton distributions at Q = m-c = 1.3 GeV; (a) 
u-quark; and (b) d-quark. The axes are scaled to highlight the valence components of these distributions. 




FIG. 4: Same as Fig. |3| except the axes are scaled to highlight the sea components of the PDFs: (a) u-quark; and 
(b) d-quark. 



distributions for the HQ scheme since these distributions are compensating for the suppression of charm in 
this scheme; this suppression is absent in the zero-mass MS scheme. In the region of ultra-small x, there 
is no data constraining these distributions; we show this effect graphically in the next Section where we 
consider the PDF uncertainty band of the the HQ and MS distributions. 

Fig. El shows the comparison of the gluon distribution at Qq. Here the difference between C5HQ and 
the CTEQ6 generation of gluon distributions is pronounced. As discussed in Ref. 0, the change in this 
least-well-determined parton distribution is due to the recent precision DIS data (most influential in the 
small X region) in conjunction with the greatly improved inclusive jet data from the Tevatron (critical for 
the medium to large x regions). The differences between C6HQ and C6M gluons at large x are due to a 
combination of scheme-dependence, and the inherent uncertainty range of the current analysis. [l| 

Fig. El shows the comparison of the strange distributions at the same Qo- The noticeable difference 
between the C5HQ curve and the others, in this case, is largely the result of different theoretical inputs: 
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CTEQ6HQ Gluon Distribution compared to other CTEQ sets 




FIG. 5: Comparison of the gluon distributions at Qo~ 1-3 GeV. 



namely, the k parameter (cf., Eq. ISJ. k determines the ratio of strange to non-strange sea quarks at the 
initial scale Qo- This k factor, known only approximately, was chosen to be 1/2 both in the CTEQ5 and 
CTEQ6 analyses, but for slightly different values of Qq—I.O GeV for C5HQ and 1.3 GcV for CTEQ6 sets. 




FIG. 6; Comparison of the strange distributions at Qo= 1-3 GeV. 
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FIG. 7: Comparison of the charm and gluon distribution at — 10 GeV^. 



B. Charm and gluon distributions at — 10 GeV^ 

To compare the difFerences of the charm distributions, we need to move some distance above the 
charm mass scale, since all current sets assume a zero charm distribution at /i = m^. Fig. [7|i makes this 
comparison at = 10 GeV^. We see a substantial difference between the C5HQ and the two CTEQ6 
charm distributions. Unlike the case of the strange quark, this difference is physical; it is mainly a reflection 
of the difference in the gluon distributions (shown already for a lower Q value in Fig. [SI because charm is 
radiatively generated from the gluon. To confirm this hypothesis. Fig. [71d shows the gluon distributions at 
the same = 10 GeV^; the similarity is clear. 

C. Q-dependence of heavy quark mass effects 

In this subsection, we compare the Q dependence of the C6M and C6HQ parton distributions at fixed 
values of x in order to see how the differences between PDFs in the zero-mass and general-mass schemes 
vary with increasing /x-scale. We expect these differences (relative to the general PDF uncertainties) will 
decrease with increasing Q. The purpose of this study is two-fold. First, we check the self-consistency of 
the analysis, which is based on the expectation that power suppressed mass terms vanish asymptotically as 
Q ^ oo. Secondly, we determine the specific scales where the mass-effects become insignificant, in practice. 
The latter question is relevant for high-energy applications of PDFs derived from lower-energy data.^'' It can 
also be seen as a general illustrative example demonstrating how the relative importance of power suppressed 
terms decrease in the high region where the conventional zero-mass PQCD results becomes dominant. 

To examine the effects due to heavy flavor masses in the theoretical formalism, we cannot directly 
compare physical quantities that are being fit — they are the same (within the errors of the global fit) by 
construction. On the other hand, parton distribution function of individual flavors are not a good gauge 
either, because they are not physical — they are scheme dependent. Thus, for this comparison, we compromise 
and examine the quark singlet combination (S = + ?)) ^.nd the gluon distribution. 



It also serves to validate the self-consistency of analyses where threshold effects have been neglected (such as for collider 
processes). 
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FIG. 8: Ratio of GM (C6HQ) to ZM (C6M) parton distributions as a function of Q [GcV] along with the un- 
certainty band of the latter. The offset on the y-axis is arbitrary; i.e. the horizontal lines all correspond to 
'£{x,Q)^'^/T,{x,Qf^ = 1. Plotted are ratios for x = {10~^ 10""*, 10"^, 10"^, 0.1, 0.5} where x is increasing 
from the lowest to the uppermost curve. We display: (a) the gluon; and (b) the singlet quark E = X]q(9 + 



Figs. |H^,b show, as solid lines, the C6HQ gluon and quark singlet distributions normalized to that of 
C6M. The shaded bands represent the ranges of uncertainty of the relevant quantities due to experimental 
sources, as estimated in the CTEQ6 global analysis. 1] We observe the following pattern: for most of 
the parameter space, the deviation of the C6HQ from C6M falls within the estimated uncertainty range. 
In particular, for the gluon distribution this is true for all values of x. For the quark distribution, the 
HQ fit overshoots the zero-mass fit significantly (i.e., beyond the uncertainty band) at moderately low x 
~ [10~^, 10"'^]. At very low x ~ 10~^, the differences fall within the uncertainty bands again (primarily due 
to the increasing uncertainty). Although the HQ quark PDFs at x ^ [10~^, 10~^] eventually evolve into the 
error band of the zero-mass fit as Q increases, one may wonder why the two fits do not merge faster. The 
GM-VFNS and ZM-VFNS differ by the manner in which they organize terms of order 0{m'^ /Q^). Thus, the 
Wilson coefficients in the two schemes will be comparable when [m? /Q'^) << 1. The boundary conditions 
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of the PDFs at Qo will also differ between the two schemes by terms of order 0{rrP' / Q^); however, since the 
PDFs evolve logarithmically in the scaling variable, while the power corrections fall off more quickly, this 
can contribute to scheme differences at scales significantly above the scale m. This effect is evident in the 
case of the singlet quark PDF at moderately low a; values. Thus, the observed behavior is quite natural. 

VI. CONCLUDING REMARKS 

The CTEQ6HQ PDFs presented here complement the previously published CTEQ6 sets by providing 
distributions which can be used in the generalized MS-bar scheme with non-zero mass partons. This analysis 
includes the complete set of NLO processes including the real and virtual quark-initiated terms. Additionally, 
the ACOT(x) scheme is used to introduce a generalized scaling variable which provides numerically stable 
results for the entire energy range — from heavy quark thresholds to the high energy limit. 

While the zero-mass parton scheme is sufficient for many purposes, the fully massive scheme can be 
important when physical quantities are sufficiently sensitive to heavy quark contributions. This is evident 
when comparing the CTEQ6HQ and CTEQ6M fits to the mis-matched sets (Table 1) where the precise DIS 
data from HERA highlights the discrepancies. 

The CTEQ6HQ fits also provide the basis for a series of further studies involving more quantitative 
analysis of strange, charm, and bottom quark distributions inside the nucleon. For example, the CTEQ6HQ 
PDFs are necessary for a consistent analysis of resummed differential distributions for heavy quark production 
such as in Ref. |37|. Using the full range of data from both the charged and neutral current processes, these 
distributions can reduce the uncertainties in the calculations; hence, they have significant implications for 
charm and bottom production, and can help resolve questions about intrinsic heavy quark constituents inside 
the proton, the Axi^s structure function, and the extraction of sinOw. 
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